Higgs boson decays in the Complex MSSM 
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Abstract. The analysis of the Higgs search results at LEP showed that a part of the MSSM 
parameter space with non-zero complex phases could not be excluded, where the lightest neutral 
Higgs boson, hi, has a mass of only about 45 GeV and the second lightest neutral Higgs boson, 
/i2, has a sizable branching fraction into a pair of hi states. Full one- loop results for the Higgs 
cascade decay h2 — > hihi are presented and combined with two-loop Higgs propagator corrections 
taken from the program FeynHiggs. Using the improved theoretical prediction to analyse the limits 
on topological cross sections obtained at LEP, the existence of an unexcluded region at low Higgs 
mEiss is confirmed. The effect of the genuine vertex corrections on the size and shape of this region 
is discussed. 

PACS. 12.60.Jv Supersymmetric models - 14.80.Cp Non-standard-model Higgs bosons 



1 Introduction 

In its general form, the MSSM contains complex phases 
which drive CT^-violation. The results of the Higgs- 
boson searches at LEP were analysed in a series of 
MSSM benchmark scenarios [1,2], including one ex- 
plicitly chosen to investigate the effects of these CV- 
violating phases, called the CPX scenario [2]. This sce- 
nario produced an unexcluded region for light hi and 
relatively small values of tan /? (the ratio of the vacuum 
expectation values of the two Higgs doublets), so that 
no firm lower bound on the mass of the lightest Higgs 
boson of the MSSM could be set [3]. Throughout much 
of this region of CPX parameter space, the process 
^2 hihi dominates the decay width. Therefore, 
in order to reliably determine which parameter regions 
of the MSSM are unexcluded by the Higgs searches 
so far and which regions will be accessible by Higgs 
searches in the future, precise predictions for this pro- 
cess are indispensable. 

However, the genuine vertex corrections (as op- 
posed to propagator corrections) to the process /12 ^ 
hihi in the MSSM with complex parameters have so 
far been unavailable within the Feynman-diagrammatic 
(FD) approach. These vertex corrections are expected 
to be sizable as they contain terms depending on the 
fourth power of the top mass. We have obtained com- 
plete one-loop results within the FD approach for the 
decays ha — > hbhc, including all genuine one- loop ver- 
tex corrections, in the MSSM with complex phases 
(cMSSM) [4]. We have furthermore calculated com- 
plete one-loop results for the decays of neutral Higgs 
bosons into SM fermions in the cMSSM [4]. These new 
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one-loop results are combined with existing higher- 
order corrections in the FD approach and will be in- 
cluded in the public code FeynHiggs [5-8] . As an appli- 
cation of our improved result, we analyse the coverage 
of the LEP exclusion limits on the topological cross 
sections for the various Higgs-boson production and 
decay cross sections within the Mh^^-tanP parameter 
plane. 



2 Calculation of Higgs-boson cascade 
decays and decays into SM fermions 

We calculated the full IPI (one-particle irreducible) 
one-loop vertex contributions to the processes ha 
hbhc, taking into account all sectors of the MSSM 
and the complete dependence on the complex phases 
of the supersymmetric parameters. These calculations 
(which make use of the programs FeynArts [9] and 
FormCalc [10]) have been described in more detail in 
Rcf. [4]. 

We have derived also complete one-loop results for 
the processes ha ff (including SM-type QED and, 
where appropriate, QCD corrections) for arbitrary val- 
ues of all complex phases of the supersymmetric pa- 
rameters, as these decays are important over large 
parts of the cMSSM parameter space. The partial de- 
cay widths for the other Higgs-boson decay modes have 
been taken from the program FeynHiggs. 

In our numerical analysis below, we will compare 
our full result with the contribution from just the t, t 
sector in the approximation where the gauge couplings 
are neglected and the diagrams are evaluated at zero 
external momenta, which we call the 'Yukawa Approx- 
imation'. 
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For the renormalisation of the Higgs fields it is con- 
venient to use a DR scheme as in Ref. [5]. Therefore, 
it is necessary to include finite wave-function normal- 
isation factors (Z-factors) to ensure that S-matrix el- 
ements involving external Higgs fields have the cor- 
rect on-shell properties. The form of these Z-factors is 
given in Ref. [4]. In the calculation of the Z-factors, we 
incorporate higher-order contributions obtained from 
the neutral Higgs-boson self-energies of the program 
FeynHiggs. In the analysis below, we will restrict to 
those higher-order contributions for which the phase 
dependence at the two-loop is explicitly known [6], i.e. 
we do not include contributions that have been extrap- 
olated from results obtained for the real MSSM. 



3 Implementation of exclusion bounds 
from the LEP Higgs searches 

In our numerical analysis below we analyse the impact 
of our new result on the LEP exclusion bounds in the 
cMSSM. This is done by comparing the cMSSM pre- 
dictions with the topological cross section limits given 
in Refs. [3, 11]. This procedure is described in more de- 
tail in Ref. [4] and Ref. [12]. It involves checking which 
channel is predicted to have the highest statistical sen- 
sitivity and then comparing the theoretical prediction 
in that channel only with the topological cross section 
limit determined at LEP. This ensures a correct sta- 
tistical interpretation of the overall exclusion limit at 
the 95% C.L.. However, it should be noted that, since 
only one channel can be used at a time, this method 
leads to poorer coverage in areas where two or more 
channels have a similar statistical sensitivity than, for 
example, the dedicated analyses carried out in Ref. [3]. 

4 Numerical results 

In our numerical analysis we use the parameter val- 
ues of the CPX benchmark scenario [2], adapted to 
the latest experimental central value of the top-quark 
mass [13] and using an on-shell value for the abso- 
lute value of the trilinear couplings At and Ai, that is 
somewhat shifted compared to the DR value specified 
in Ref. [2]. Specifically, if not indicated differently, we 
use the parameters 

A^susY ^ 500 GeV, \At\ = \Ab\ = 900 GeV, 

fi = 2000 GeV, rug = [il/g] ^ 1000 GeV, 

mt = 170.9 GeV, Ah = 200 GeV. (1) 

The lowest-order Higgs sector parameters tan /3 and 
Mh± are varied, with A/^± < 1000. The complex 
phases of the trilinear couplings At, Ah and the gluino 
mass parameter are set to 

TT TT TT 

VAt = 2 , ^A, = 2 , ^~g^ 2 • (2) 

In Eq. ([T]) MgusY denotes the diagonal soft SUSY- 
breaking parameters in the sfermion mass matrices, 
which are chosen to be equal to each other. 
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Fig. 1. r{}i2 h\h\) as function of tan/? {Mu± is ad- 
justed to ensure Mh^ = 30 GeV), for various approxima- 
tions. 
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Fig. 2. The branching ratio BR(/i2 — > h\li\) in the M^^- 
tan P plane of the CPX scenario. 

4.1 Results for the h2 hihi decay width 

Fig. [T] shows the relative effect of different contribu- 
tions to the /i2 hihi decay width in an area of the 
cMSSM parameter space that is particularly relevant 
to the unexcluded region in the LEP Higgs searches. 
The decay width is plotted against tan/3 whilst ad- 
justing Mh± such that Mh^ = 30 GeV. All the re- 
sults plotted include the higher-order corrected wave- 
function normalisation factors as described in Sect. [3) 
The results differ only in the genuine contributions 
to the h2hihi vertex. One can see that the full re- 
sult (denoted as 'Full') differs drastically from the case 
where only wave-function normalisation factors but no 
genuine one-loop vertex contributions are taken into 
account ('Tree'). The Yukawa approximation agrees 
much better with the full result and explains the qual- 
itative behaviour of the full result — in particular, the 
region where r{h2 — > hihi) « for tan/3 « 4.3 is 
due to the fact that the Yukawa vertex corrections to 
the matrix element change sign when tan/3 is varied. 
Using the full contribution from the t, i, 6, h sector ('t, 
st, b, sb') or using all three generations of fermions 
and sfermions ('f, sf) yields a prediction that deviates 
from the full result by up to 15%. 
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Fig. 3. The channels predicted to have the highest sta- 
tistical sensitivity for setting an exclusion limit using the 
results from the LEP Higgs searches in the CPX scenario. 
The colour codings are: red = (hiZ) — > (bbZ), blue — 
{hiZ) -» {bbZ), white = {h2Z) -* (hihiZ) {bbbbZ), 
cyan = (/12/ii) ~^ (bbbb), yellow — (h2hi) —* (hihihi) 
(bbbbbb), black = other channels. 
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Fig. 4. The parameter region in the CPX scenario ex- 
cluded at the 95% C.L. by the topological cross section 
limits obtained at LEP [3]. The colour codings are: green 
— LEP excluded, white — LEP allowed, grey — theoreti- 
cally inaccessible. 



In Figs. [21 [3] and [4] the Mh^-tanP parameter space 
of the CPX scenario is analysed. Fig. [2]shows the branch- 
ing ratio of the Higgs cascade decay, BR(/i2 — > hihi). 
Over a large part of the parameter space where the 
decay /12 — > hihi is kinematically possible it is actu- 
ally the dominant decay channel. The branching ratio 
is particularly large for low and moderate values of 
tan/3. In the region where tan/? k, 4-5 the Yukawa 
contribution to the matrix element for the /12 hihi 
decay changes sign and causes a sharp drop in the 
/i2 hihi branching ratio, as was already observed 
in Fig. [TJ A similar behaviour occurs also in the region 
of tan P w 35. 

Fig. [3] indicates which channel has the highest sta- 
tistical sensitivity and therefore which channel will be 
used to set an exclusion limit in different regions of the 



parameter space. As explained in Sect.[4l this infor- 
mation is needed for a statistical interpretation of the 
topological cross section limits obtained at LEP. One 
can see in the figure that the channels e'^e^ — > —>■ 
hihiZ bbbbZ and e~^e~ /12/11 — > hihihi 
bbbbbb have the highest search sensitivity in a region 
with small M^^ and small and moderate values of tan {3 
(the region at tan/3 » 4 where the channel e+e^ 
(/i2^) — > (bbZ) has the highest search sensitivity cor- 
responds to the drop in the /12 ~* hihi branching ra- 
tio observed in Fig. [2]). For small Mji-^ and somewhat 
higher tan/3 the channel e+e^ (^2^1) — > (bbbb) has 
the highest search sensitivity. It should be noted that 
all channels involving the decay of the /12 boson in 
the region of small are strongly influenced by the 
r{h2 —> hihi) decay width, either directly in the case 
of the channels involving the Higgs cascade decay, or 
indirectly through the branching ratio of the /i2. The 
parameter region where /^(/i2 ^ ^1^1) is important 
coincides with the region of the CPX scenario that 
could not be excluded at the 95% C.L. in the analysis 
of the four LEP collaborations [3]. 

In Fig. m we have compared our new theoretical 
predictions with the topological cross section limits 
obtained at LEP for the channels in Fig. [3l We find 
an unexcluded region at Mh^ » 45 GeV and moderate 
tan/3 where channels involving the decay /12 hihi 
play an important role. Thus, our analysis, based on 
the most up-to-date theory prediction for the /12 — > 
hihi channel, confirms the 'hole' in the LEP cover- 
age observed in Ref. [3] (see in particular Fig. 19 of 
Ref. [3]). 

It should be noted, on the other hand, that the re- 
sults shown in Fig. [4] differ in several respects from the 
results presented in Ref. [3]. As discussed above, our 
analysis has less statistical sensitivity near borders be- 
tween areas where different search topologies are pre- 
dicted to have the highest exclusion power than the 
benchmark studies of Ref. [3] . A further difference (be- 
sides the improved theoretical prediction) is the input 
value of the top-quark mass. While we are using the 
latest experimental central value of = 170.9 [13], 
most of the analysis of Ref. [3] was done for nit = 
174.3. We have exphcitly checked that (as expected) 
the unexcluded region is significantly increased if we 
use nit = 174.3 instead. It should also be noted that 
the shape and size of the unexcluded region depends 
sensitively on the value of At used. 

In Fig. Owe focus on the parameter region at ~ 
45 GeV and moderate tan f3. While plot (a) shows the 
full result (i.e., it is a detailed view of Fig. |4]), in plot 
(b) the genuine /i2 ^ hihi vertex corrections arc ap- 
proximated by the contributions from fermions and 
sfermions only, and in plot (c) the Yukawa Approxi- 
mation has been used for the genuine /12 h\hi ver- 
tex corrections. In all three plots the wave function 
normalisation factors and all other decay widths are 
the same. While all three graphs show the unexcluded 
region to be in a similar position in parameter space, 
the shape of this unexcluded region changes quite sig- 
nificantly. 
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Fig. 5. A closer look at the LEP coverage for M^^ « 45 GeV and moderate tan/3. Plot (a) shows the full result (detailed 
view of plot (b) of Fig. |3} . Plot (b) shows the result for the case where only contributions from SM fermions and their 
superpartners are taken into account in the genuine vertex corrections. Plot (c) shows the result where the Yukawa 
Approximation has been used for the genuine vertex corrections. The colour codings are: green = LEP excluded, white 
= LEP aUowed. 



5 Conclusions 

We have shown that our new result for the decay width 
ft,2 — > hihi, which includes the complete one- loop ver- 
tex corrections, drastically improves on results obtained 
using propagator corrections alone. We have combined 
this with new results for neutral Higgs decays into two 
fermions in order to re-analyse the CPX parameter 
space using the limits on topological cross sections ob- 
tained from the LEP Higgs searches. We find that over 
a large part of the parameter space where the decay 
ft,2 — > hihi is kinematically possible, it is the dominant 
decay channel. We confirm that the parameter space 
at low Higgs mass can not be completely excluded by 
the LEP Higgs searches. We furthermore show that the 
size and shape of the unexcluded region is significantly 
modified if approximations to the full result are used, 
emphasizing the relevance of the full one-loop result 
for the h2 hihi decay width in the complex MSSM. 
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